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We introduce and develop a novel approach to extend the ordinary two-flavor neutrino oscillation 
formalism in matter using a non-Hermitian PT symmetric effective Hamiltonian. The condition of 
PT symmetry is weaker and less mathematical than that of hermicity, but more physical, and such 
an extension of the formalism can give rise to sub-leading effects in neutrino flavor transitions similar 
to the effects by so-called non-standard neutrino interactions. We derive the necessary conditions for 
the spectrum of the effective Hamiltonian to be real as well as the mappings between the fundamental 
and effective parameters. We find that the real spectrum of the effective Hamiltonian will depend 
on all new fundamental parameters introduced in the non-Hermitian PT symmetric extension of 
the usual neutrino oscillation formalism and that either i) the spectrum is exact and the effective 
leptonic mixing must always be maximal or ii) the spectrum is approximate and all new fundamental 
parameters must be small. 
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In 1998, the Super-Kamiokande collaboration [l| re¬ 
ported results that gave a clear indication of neutrino os¬ 
cillations —a quantum mechanical effect over distances 
of thousands of kilometers, and thus the first solid ev¬ 
idence for physics beyond the Standard Model of par¬ 
ticle physics. Now, it has been firmly established that 
neutrino oscillations is the best and leading description 
of neutrino flavor transitions and Takaaki Kajita and 
Arthur B. McDonald have been awarded the Nobel Prize 
in Physics 2015 “for the discovery of neutrino oscillations 
... ”. However, other mechanisms could be responsible 
for such transitions on a sub-leading level. Therefore, 
other “new physics” effects have been proposed, e.g. so- 
called non-standard neutrino interactions (NSIs). See, 
e.g., ref. Q and references therein. 

Nevertheless, in the same year as the Super- 
Kamiokande collaboration reported its results, Bender 
and Boettcher Q presented (based on a conjecture by 
Bessis) PT symmetric quantum mechanics, which is de¬ 
scribed by non-Hermitian Hamiltonians. In ordinary 
quantum mechanics, a Hamiltonian operator H should be 
Hermitian (ijl = H) m order to have real (and measur¬ 
able) energy eigenvalues. However, using PT symmetric 
quantum mechanics, the requirement of hermicity can be 
replaced by a weaker, but more physical, requirement of 
space-time reflection symmetry (so-called PT symmetry, 
i.e. [PT,i?] = 0) without losing any important physical 
aspects of ordinary quantum mechanics y . Note that P 
is parity (space reflection) and T is time reversal. Now, 
PT is an antiunitary operator . The definition for any 
such operator A is that (A(()|A'0) = {(j)\ijj)*, where \ip) 
and {(j)) are arbitrary states. Any antiunitary operator 
can be written in the form A = UK, where U is a, uni¬ 
tary operator {e.g., P) and A' is a complex conjugation 
{e.g., T in position representation). Consider Hamilto¬ 
nians with antiunitary symmetry, i.e. [A, H] = 0, where 
^ 2 k = X (fc odd). Note that PT corresponds to fe = 1. It 


has be shown that for any such A, it is possible to con¬ 
struct a basis in which the matrix elements of K are real, 
which means that the characteristic equation giving the 
eigenvalues of II is real Q . This result is a generalization 
of applications in quantum chaology Q, which, in turn, 
was a generalization of arguments from nuclear physics 

i. 

In refs. , it has been shown that two Dirac 

fermions coupled by a non-Hermitian PT symmetric mass 
matrix describes a single 8-dimensional relativistic par¬ 
ticle of a fundamentally new type that can appear as 
two massless particles, despite the mass matrix being 
non-zero. It is then evident that a non-Hermitian PT 
symmetric mass matrix might lead to a different neu¬ 
trino oscillation phenomenology, since such a mass ma¬ 
trix could describe neutrinos that oscillate between two 
flavors (given by the two massless particles) but prop¬ 
agate masslessly. Now, evidence of neutrino oscilla¬ 
tions would normally mean that neutrinos are massive. 
However, in PT symmetric quantum mechanics, the 8- 
dimensional solution suggests that this conclusion might 
not need to be drawn. Later, in ref. till, it has been 
discussed that interactions between ordinary light neu¬ 
trino states and new very light many-particle states could 
lead to new physics described by non-Hermitian Hamil¬ 
tonians. The idea has been developed for neutral kaons 
and can be readily adapted to neutrinos, which could 
interact with and decay into such generic light states. 
For neutrinos, there are no constraints from CPT in¬ 
variance, whereas for neutral kaons, this is not the case. 
Recently, in ref. 12|, a non-Hermitian Yukawa theory 
has been studied that could give an explanation for the 
smallness of the masses of light neutrinos. Finally, con¬ 
necting PT symmetric quantum mechanics to neutrino 
oscillations. Bender et al. have experimentally studied a 
simple mechanical system consisting of two coupled pen- 
dula, which is the classical analog of the phenomenon of 
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neutrino oscillations. They have observed a phase tran¬ 
sition that separates the unbroken and broken phases of 
the PT symmetry of the quantum-mechanical Hamilto¬ 
nian [l^ . Therefore, this simple experiment and its result 
provide an intuitive motivation to study non-Hermitian 
PT symmetric Hamiltonians for neutrino flavor transi¬ 
tions in matter. 

In this Letter, we will introduce and investigate two- 
flavor neutrino oscillations in matter (of constant den¬ 
sity) based on non-Hermitian, but PT symmetric, Hamil¬ 
tonians that could give rise to sub-leading effects in neu¬ 
trino flavor transitions. Based on the physical motiva¬ 
tion of earlier results [Ml, and especially the result of 
ref. 0 , such an extension of ordinary neutrino oscilla¬ 
tions will therefore open up the possibility for new physics 
beyond the Standard Model. In general, non-Hermitian 
extensions relax assumptions about constraints made on 
Hermitian models (such as the requirement of hermicity) 
rather than adding new dynamical degrees of freedom. 
Thus, they constitute a different paradigm of new physics 
compared to traditional model building. 

Assuming two lepton flavors (for generality, they will 
be denoted a and /3), the time evolution of the neutrino 
vector of state v = {ya describing neutrino oscilla¬ 

tions is given by a Schrodinger-like equation with a time- 
independent Hermitian Hamiltonian Hose, namely 


dt 


I 
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2 
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- Hosc^(0 , 


( 1 ) 


where E is the neutrino energy and A = 2y/2EGFNe 
is the effective matter potential induced by ordi nary 
charged-current weak interactions with electrons [IJ, ll5| . 
Here, mi and m 2 are the definite masses of the neu¬ 
trino mass eigenstates, vi and 1 ^ 2 , respectively, that 
are related to the weak interaction eigenstates, Va and 
Fp, through the leptonic mixing matrix O such that 

(t'o = O 1 ^ 2 )"^, Gf is the Fermi coupling con¬ 

stant, and Nf. is the electron density of matter along the 
neutrino trajectory. Furthermore, the two-flavor leptonic 
mixing matrix (parametrized by the leptonic mixing an¬ 
gle 9, i.e. one real parameter) can be written as 


O = 


—s c 


€ SO{2), 


( 2 ) 


oscillations probability formulas are given by (see, e. 5 ., 
ref. 161) 


P{va vp-,L) = sin^(20) sin^ 

P{Va —^ r'ai A) = 1 — P{Va Vp] L) 

= 1 — P{vp —>■ r'a! L) = P{i'p —>■ yp] L) , (4) 


where L is the (propagation) path length of the neu¬ 
trinos and Am^ = is the mass-squared dif¬ 

ference between the masses of the two neutrino mass 
eigenstates. Note that 6 corresponds to the amplitude 
of the oscillations, whereas Am^ corresponds to the fre¬ 
quency. Indeed, in matter of constant density {i.e. as¬ 
suming Ne = const, i 0) as well as in the case of NSIs, 
the probability formulas are obtained by replacing the 
vacuum parameters with effective matter or NSI para¬ 
meters in eqs. m and dl]). (Note that if Ng = const., 
A o<: E, and hence, A is only constant for a fixed E.) 
Thus, we have the transition probability 

/ /\m'^ T \ 

P{yo, ^ yp-, L) = sm\2e') sin^ , (5) 


where 9' is the effective mixing angle and Am'^ is the 
effective mass-squared difference. Therefore, there are 
non-trivial mappings between the vacuum parameters 9 
and Ato^ and the corresponding effective (matter or NSI) 
parameters 9' and Am'^. The key point will be the di- 
agonalization of a given effective Hamiltonian Heff and 
obtaining the explicit relations of the effective parame¬ 
ters in terms of the fundamental (vacuum) parameters 
0 . 

In the case of matter of constant density, the mappings 
are given by 


sin2(26»') = 


sin^( 20 ) 


sin^(2d) -I- (cos(20) — A/Am^) 


( 6 ) 


Am'^ = \J (Am^)^ sin^( 20 ) -|- (Am^ cos( 2 d) — A)^ 


( 7 ) 


Note that the effective matter mixing angle in eq. ([HI) ob¬ 
tains its maximal value sin^( 20 ') = 1 when the condition 

A = Ami cos( 2 d) ( 8 ) 


where c = cos 9 and s = sin0. Using the formalism 
of neutrino oscillations, the quantum mechanical transi¬ 
tion probability amplitudes are given as overlaps of differ¬ 
ent neutrino states, and eventually, neutrino oscillation 
probabilities are defined as squared absolute values of 
the amplitudes. Thus, neutrino flavor transitions occur 
during the evolution of neutrinos. In vacuum (i.e. as¬ 
suming Ng = 0 (or A = 0 )), using conservation of proba¬ 
bility (i.e. unitarity), the well-known two-flavor neutrino 


is fulfilled. This condition is the famous Mikheyev- 
Smirnov-Wolfenstein (MSW) resonance condition [l^, 


18|. If the MSW condition is satisfied, then the effective 


matter mixing angle 9' is maximal {9' = 45°) indepen¬ 
dently of the value of the fundamental vacuum mixing 
angle 9. Thus, the transition probability in eq. can 
be very large even if 9 is very small. In the vacuum 
limit, i.e. the limit of vanishing A, the effective parame¬ 
ters are replaced by the fundamental vacuum parameters 
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as well as eqs. m and (lU are identical to each other and 
the mappings ([6]) and 0 are reduced to 6' = 0 and 
Am'^ = I Ato^ I, respectively. 

In order to study the features of adding a non- 
Hermitian, but PT symmetric, Hamiltonian to the or¬ 
dinary formalism of neutrino oscillations in matter, we 
consider for illustration the following general complex 
Hamiltonian UM- 





(9) 


where p, a, and p are three real parameters. The physical 
meaning of p, cr, and p is similar to other phenomenologi¬ 
cal parameters for “new physics” {e.g., NSI parameters). 
The Hamiltonian H' is not Hermitian, but it is PT sym¬ 
metric. However, note that if (p = 0,7r, then H' be¬ 
comes Hermitian. Using the characteristic equation of 
this Hamiltonian, i.e. <let{2EH' — AI 2 ) = 0, we find the 
two eigenvalues A = Ap cos p ± H\/cr^ — p^ sin^ p. Thus, 
there are two parametric regions. If sin^ p, the 

two eigenvalues form a complex conjugate pair and the 
PT symmetry is broken, whereas if > p^ sin^ p, the 
two eigenvalues are real and the PT symmetry is unbro¬ 
ken. At the point between the two regions, i.e. = 
p^ sin^ p, there is only one eigenvalue and it is called an 
exceptional point. The interesting region to consider is 
the unbroken region with conservation of PT symmetry 
and two real eigenvalues. Note that complex Hamilto¬ 
nians may be either Hermitian or PT symmetric, but 
not both, whereas real symmetric Hamiltonians may be 
both Hermitian and PT symmetric Q. Another way 
of introducing a complex non-Hermitian Hamiltonian in 
a different context has been performed for high-energy 
two-flavor neutrino oscillations in matter with absorption 
leading to complex indices of refraction [l9| . 

Now, constructing the effective Hamiltonian Hgg by 


adding the non-Hermitian PT symmetric Hamiltonian 
H' to the ordinary neutrino oscillation Hamiltonian in 
matter Hose and following refs. [3 [23 , we obtain the 
total Hamiltonian, which is complex symmetric, as 


Hefi — Hqsc 
_ 1 


H' 


O 


0 m 2 




pe 


+ A 


2E^ 


a 

-lip 

0 

0 m" 


pe 


9 

mf 




A 0 
0 0 
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where C/ is the effective leptonic mixing matrix and mi 
and m 2 are the two effective neutrino masses, which are 
assumed to be non-degenerate, otherwise there would not 
be any oscillations. Since H^s is a 2 x 2 complex sym¬ 
metric matrix, U must be a 2 x 2 unitary matrix such 
that 


ii = 



se'^2 

ce-‘‘^i 


e su{2 ), 


( 11 ) 


where c = cos 0, s = sin 9, 9 being the effective mixing an¬ 
gle, and Pi and p 2 are two phases, i. e. in total three real 
parameters. Note that the total Hamiltonian can be ex¬ 
pressed in three different useful and relevant bases of the 
same Hilbert space, namely the flavor, mass, and effec¬ 
tive bases. Thus, we have = UHeU'^ = OHmO~^, 
where and Hm are the total Hamiltonian in effective 
basis and mass basis, respectively. In the case of two lep¬ 
ton flavors, O is orthogonal, and therefore, 0~^ = . 

Thus, diagonalizing the effective total Hamiltonian HeS 
in eq. OT, we have the following six relations: 


c^mi -I- s^mf -I- A -I- Apeosp = cos(2pi)mi -I- cos(2p2)mf , (12) 

Apsinp = sin(2pi)m^ -|- sin(2p2)m2 , (13) 

scAm^ -I- Aa = sccos(pi — p 2 )Am^ , (14) 

0 = scsin(pi — p 2 )Am^ , (15) 

(?m\ -\- Apeosp = cos(2p2)mi -I- (? cos(2pi)mf , (16) 

Apsinp = sin(2p2)m2 + c^sin(2pi)m2 , (17) 


where Am^ = m^ — mf. Note that eq. (IT5|) immediately 
implies that pi = p 2 . Solving eqs. (inii-diTl) for the 
effective amplitude and frequency, we find the possible 


mappings 


sin^(2d) 


(Am^ sin(20) -|- 2Acr)^ 

(Am^ sin(2d) -|- 2Aa)^ -\- {Am? cos(2d) — A)^ 

(18) 


An? 


\J (Am^ sin(20) -|- 2 Act)^ -|- {Am? cos(2d) — A)^ . 


(19) 
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In addition, we obtain two auxiliary mappings 


ference, respectively, are reduced to 


tan(2(^) 


~ 2 I ~ 2 

rrii + TO2 

where (p = In the vacuum limit [i.e. A —> 0), 

for eqs. (fT^ - (l^ . we have 6 = 6, Am^ = |Am^|, p = 0, 
and rhi + rh^ = \mf + mf]. In this case, note that the 
effective phase p naturally becomes equal to zero, since 
it does not have a correspondence in ordinary vacuum 
neutrino oscillations. 

Does the effective total Hamiltonian Hgg have a re¬ 
gion of unbroken PT symmetry? In order to answer this 
question, we need to consider the characteristic equation 
of this Hamiltonian, i.e. det{2EHes — AI 2 ) = 0, which 
leads to a quadratic equation 


2ApsiTLp 


-I- TO 2 + ^ + cos ip 


( 20 ) 


\J{m\ + m\ + A + 2Ap cos p)^ + dA^p"^ sin^ p , 

( 21 ) 


sin2(20) = l, (25) 

Am^ = |Am^sin(26*) -|-2 Hct| > \2Aps\np\ . (26) 

We observe that the value of the effective mixing angle 
is always maximal [6 = 45°) independently of the values 
of the three parameters p, a, and p describing H'. The 
reason is that the MSW condition has to be satisfied in 
eq. (UHl). Furthermore, we note that the effective mass- 
squared difference is always positive and only linearly 
dependent on the parameter cr. 

The effects of a non-Hermitian PT symmetric Hamilto¬ 
nian on neutrino flavor transitions must be sub-leading, 
and therefore, the parameters describing H' in eq. (|^ 
have to be small. Thus, we series expand the two aux¬ 
iliary mappings o and (1^ up to second order in the 
small parameter p (note that the parameter is a phase): 


— {mf + 1711 + A + 2Ap cos p)X 

+ {c^ml + s^ml + A){s^ml + c^m^) 

+ A^p^ — {scAm? + AaY 
+ {s^m\ + c^ml)Ape'‘^ 

+ {c^m\ + + A)Ape~"^ = 0 , ( 22 ) 


where the last two terms on the left-hand side of this 
equation contain imaginary parts. The condition for 
the equation to have real roots is that the imaginary 
parts equal zero, see, e.g., ref. [21], pp. 91-92. It 
turns out that the imaginary parts can be reduced to 
S = Apsmp{Ani^ cos(20) — A), which means that 3 = 0 
if the MSW condition in eq. m is fulfilled. Of course, 
3 = 0 when either p = 0 or p = 0,Tr, but these are triv¬ 
ial cases and therefore not interesting. (Later, we will 
investigate another non-trivial, but approximate, case.) 
Thus, the MSW condition is a necessary constraint for 
having two potential real exact eigenvalues of iLeff- 
Inserting the MSW condtion ([U into eq. (1221) . we ob¬ 
tain 


\ 2 2 I 2 2 I ,1 

A = s nil + c 1712 + Ap cos p 


p = 


Hsin(/5 




■A 


p 0{p^), 


-I- m2 = \mi -I- m2 -I- H| 

-I- 2 sgn(mi -I- m 2 -I- A)Ap cos p + 0{p‘^), 


(27) 

(28) 


which are both directly independent of the parameter a. 
Indeed, we observe that the effective phase in eq. (EZl) has 
a sinusoidal dependence on p for small p. In addition, 
using eqs. (l26ll and (l2^ . we obtain the series expansions 
of the two real eigenvalues of H^s as 

ml = 2 ( 1 ™? ^^2 + ^1 ~ I Am^ sin(20) -|- 2 H(t|) 

sgn{ni\ni 2 +A) Ap cos p + 0{p'^), (29) 

™2 = 2 ^^2 + ^1 + I Am^ sin(20) -|- 2 H(t|) 

+ sgn{m'l-\-ni\ +A) Ap cos p + 0{p^). (30) 


Now, combining the MSW condition ([8|) and the PT 
inequality (l24)) . we obtain a single inequality among 
the fundamental parameters in order to have a non- 
Hermitian PT symmetric effective total Hamiltonian 
with a real exact spectrum, namely 


± [scAw? + AcrY — A^p'^ sin^ p . (23) 


A^ [4(psini,£5 — cr)^ + 1 ] < (Am^)^ , (31) 


In order for eq. (l23|) to have two real eigenvalues, the 
quadratic form {scAm? + Aa)^ — A^p^ sin^ p must be 
positive definite, which consecutively means that the fol¬ 
lowing condition needs to be fulfilled: 


which holds if both Am^ sin(20) -|- 2Aa and Ap s\np are 
positive or negative. In the case in which one of these 
two expressions is positive and the other one is negative 
or vice versa, the corresponding inequality is 


I Am^ sin(20)-I-2 Hct| > \2Aps\rLp\ . (24) 

Thus, iLeff has a region of unbroken PT symmetry if the 
two conditions in eqs. ([8|) and (l24)l are satisfied. In this 
case, using eqs. and (HU), the mappings for the ef¬ 
fective mixing angle and the effective mass-squared dif- 


A^ [4(psin(/5 -I- (j)^ + 1 ] < (Am^)^ . (32) 

Since Am? (or equivalently, 6 via cos(20) = AjAm? = 
2^/20An?, where Fires is the resonance energy) 
is the considered two-flavor neutrino mass-squared differ¬ 
ence (given by Nature) and A is the effective matter po- 
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FIG. 1: The spectrum of the effective mass-squareds mf and 
mi as functions of (p G [0, 2tt) with cr = 0 and mi = 10~® eV^. 
The solid curves show the case of atmospheric neutrinos 
(p = 0.1), whereas the dashed curves show the case of solar 
neutrinos (p = 1). A red curve depicts mi, whereas a green 
curve depicts mi. In addition, a blue curve displays Am^ = 
mi — mi. The fundamental neutrino oscillation parameters 
used are: Amii = 2.5 • 10“® eV^, Amii = 7.5 • 10“® eV^, 
013 = 8.5°, 012 = 33°. 


tential for the chosen experimental setup {e.g., a specific 
neutrino oscillation experiment), the model parameters 
p. O’, and if have to be restricted such that eq. m (or 
eq. (1321) ) is fulfilled. Thus, the important inequality (IMI) 
sets the bound on the new fundamental parameters that 
still allows the effective total Hamiltonian (Hoj to have 
an unbroken PT symmetry with real exact eigenvalues. 

For example, in the case of atmospheric neutrinos, we 
have the parameter values 6 ~ O 13 = 8.5°, Am? ~ 
Am §2 = 2.5 • 10“^ eV^, and ~ 11 GeV, whereas 
in the case of solar neutrinos, we have 9 ~ 612 = 33°, 
Am2 ~ Amii = 7.5 • eV^ and £1®, ~ 0.13 GeV 
HEi. Thus, inserting the two sets of parameter values 


into eq. m, we find the following estimates of the upper 
bounds on the model parameters: psintp — cr < 0.15 (for 
atmospheric neutrinos) and psinp — a ^ 1.1 (for solar 
neutrinos). 

In fig.m using eqs. m and (1^ as well as the fact that 
mi + mi + H = 2(mi + c^Am^), we present an illustration 
of the effective mass-squareds and mi, which consti¬ 
tute the spectrum of HeS, for the cases of atmospheric 
and solar neutrino oscillations. For completeness, we also 
display for the two cases the effective mass-squared dif¬ 
ference given in eq. (1211), which is independent of p, and 
are therefore shown as straight lines in the figure. 

We have shown that eqs. (I23l) - (l26l) and (l2^ - (l32l) are 
only valid at the MSW resonance described by the MSW 
condition (H)). However, note that eqs. (|27|) and (l2^ do 
not assume this condition and are therefore valid away 
from this resonance. What happens for an experimen¬ 


tal setup that does not fulfill eq. m or (|32|) ? Which 
Hamiltonian should be used away from (or at least not 
close to) the MSW resonance? In fact, the same Hamil¬ 
tonian, i.e. eq. m, should be used everywhere, but, as 
mentioned above, the effects of a non-Hermitian PT sym¬ 
metric Hamiltonian on neutrino flavor transitions must 
be sub-leading, and therefore, the new fundamental pa¬ 
rameters must be small. Assuming the three fundamen¬ 
tal parameters p, a, and ip to be small (i.e. p,a,<p 1 ) 
and keeping only terms up to first order in perturbation 
theory, the imaginary parts in eq. ( 1221 ) also equal zero, 
since 

3 = Ap sin ip{Am^ cos(26() — A) 

= ApLp{Am? cos(20) — A) 0{pp^) ~ 0, (33) 

where the factor pp is second order in the small para¬ 
meters. This is the other non-trivial case of obtaining 
a region of unbroken PT symmetry with two real eigen¬ 
values of i7eff- Therefore, in this case, eq. (1221) reduces 
to 

— {m^ m 2 A 2Ap)A 
-I- {c^m\ + a^m\ + A)(s^m^ -|- c^m\) 

— {scAm^){scAm? + 2Aa) ~ 0 . (34) 

Note that this case is only approximate, whereas the 
other case described by eq. m or (1^ is exact. How¬ 
ever, instead of solving eq. (I34|) directly, we will perform 
the same approximations in eqs. (IH-dSIl), which will 
lead to the desired results in a less tedious way. Now, 
using eq. dini) to observe that ~ 0 , since psin:/? = 
pp + 0 {pp^) ~ 0 , and series expanding eqs. (IT 2 ]l and 
(HH) in the small parameters p and p, we obtain 

(?m\ s^ml -\- A + Ap (?m\ -\- s^ml , (35) 

s mi + c m 2 + Ap ~ s mi + c m 2 , ( 35 ) 

which, combined, lead to 

ml -I- m 2 — -I- TO 2 + ^ + . (37) 

In order to find the two real approximate eigenvalues, 
i.e. the effective neutrino mass-squareds mf = ^(m^ -|- 
rhl — Am?) and m^ = 5 (mi + m^ + Air?), we need to 
compute Am^, which can be achieved by eq. (fTOl) . Series 
expanding eq. (1191) in the small parameter cr, we obtain 

An? = A"^ + (Am^)^ — 2AAm^ cos(26() 

2HAm^sin(20) 

-b , ^ ’ a + 0 {a^). 

y/ + (Am^)^ — 2 HAm^ cos(26() 

(38) 

Observing that + (Am^)^ — 2AAm? cos(26() = 

{Am? sin(20))^ -b {Am? cos(20) — A)^ and using eqs. (| 6 ]) 
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and 0 , we can therefore write eq. (1381) as 

Am^ = Am'^ + 2Asin(20')cr + O(cr^). (39) 

Thus, using eqs. (1571) and (1551) . we find the two real app¬ 
roximate eigenvalues of i7eff, 

1 2 

fh\ ~ - {m\ + m| -j- A — Am' ) + Ap — A sin(26*')CT , 

^ (40) 

7712 — + ^2 + A Am'^) + Ap + Asuy{20')<7 . 

(41) 

Furthermore, series expanding eq. m, we obtain the 
effective mixing angle, 

sin 2 ( 20 ) =sin2(26»') 

4A(2l-A7772cos(26»))2sin(26»') _ , 

+ --o- + c>(cr ), (42) 

(Atti' y 

which is consistent with eq. (l25)l if the MSW condition (|8|) 
is satisfied. Indeed, eqs. (gOl) and m are also con¬ 
sistent with eqs. ( 1291 ) and TO, since Am'^ sin( 20 ') = 
A? 77 ^ sin(20) always holds. In the vacuum limit, we of 
course have = m^, ? 7 i| = m|, and 0 = 0. Thus, when 
A for an experimental setup is not close to the MSW res¬ 
onance, we cannot use the results in eqs. (I23l) - (l26l) and 
(l29l) (l32l) . In this case, we need to use perturbation the¬ 
ory in the small parameters p, cr, and (p in order to find 
two real approximate eigenvalues of iJeff, which are given 
by eqs. TO and (TO). 

Naturally, the novelty of using a non-Hermitian PT 
symmetric effective Hamiltonian can be readily extended 
to more than two lepton flavors as well as other classes 
of the effective Hamiltonian itself. It is also beyond the 
scope of this Letter to investigate the phenomenology of 
the presented non-Hermitian PT symmetric two-flavor 
neutrino oscillation formalism, but it would certainly be 
interesting to perform such investigations for present and 
future neutrino oscillation experiments. Furthermore, it 
is obvious that the new model parameters (introduced 
by the non-Hermitian PT symmetric Hamiltonian) will 
mimick (or “fake”) effects of the old and ordinary fun¬ 
damental neutrino oscillation parameters (^ . and such 
mimicking effects could also be explored in future works. 

In summary, we have shown that a model for two-flavor 
neutrino oscillations in matter based on a non-Hermitian 
PT symmetric Hamiltonian described by three real pa¬ 
rameters p, (T, and ip has a real spectrum if either the 
condition H^[ 4 (/ 9 sini^ =F a)^ + I] < (Atti^)^ is fulfilled or 
p, cr, and p are small. In the first case, the spectrum is 
exact and the amplitude of the neutrino transition prob¬ 
ability is maximal, since the effective leptonic mixing is 
always maximal due to the fact that the ordinary MSW 
condition has to be “automatically” satisfied, whereas in 


the second case, the spectrum is approximate. In con¬ 
clusion, a PT symmetric model opens up the window to 
“new physics”, but the additional parameters need to be 
small in order for the effects to be sub-leading only, since 
neutrino oscillations are considered to be the most plau¬ 
sible description of neutrino flavor transitions. 

I am grateful to Avadh Saxena for bringing the seminal 
work on PT symmetric quantum mechanics by Bender 
and Boettcher to my attention. 
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